pi = 0
N = 100
for k in range(N):
    pi += 1 / pow(16, k)* (4 / (8*k+1) - 2/(8*k+4) - 1/(8*k+5) - 1/(8*k+6))
print("圆周率值是：{}".format(pi))